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Abstract 

Ciliates (protozoa) are ubiquitous components of plankton community and play important roles in aquatic ecosystems in 
regards of their abundance, biomass, diversity and energy turnover. Based on the stratified samples collected from the 
northern Beibu Gulf in August 201 1, species composition, abundance, biomass, diversity and spatial pattern of planktonic 
ciliates were studied. Furthermore the main environmental factors controlling ciliate communities were determined. A total 
of 101 species belonging to 44 genera and 7 orders (i.e., Oligotrichida, Haptorida, Euplotida, Sessilida, Pleurostomatida, 
Scuticociliatida and Tintinnida) were identified. The variation of ciliate communities was significant at horizontal level, but 
that was not at vertical level. Based on cluster analysis, ciliate communities were divided into three main groups. 
Redundancy analysis (RDA) revealed that Group A, existing in the waters with higher concentration of phosphorus and 
nitrogen, was dominated by Tintinnidium primitivum. Group B in the waters with lower temperature and chlorophyll-o 
concentration, was dominated by Leegaardiella ovalis. Group C, existing in the waters with higher temperature and 
chlorophyll-o concentration, was dominated by large Strombidium spp. and Mesodinium rubrum. Combining multiple 
analytic methods, our results strongly supported that phosphorus, nitrogen and chlorophyll-o were the most significant 
factors affecting the ciliate communities in the northern Beibu Gulf in summer. Concentration of phosphorus and nitrogen 
primarily influenced ciliate biomass, implying a potential impact of eutrophication on ciliate growth. The correlation with 
chlorophyll-o concentration, on one hand indicate the response of ciliates to the food availability, and on the other hand, 
the ciliates containing chloroplasts or endosymbionts may contribute greatly to the chlorophyll-o. 
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Introduction 

Ciliates are ubiquitous and significant components of micro- 
plankton communities in the ocean and play a crucial role in 
microbial food webs. They are important consumers of the pico- 
and nano- sized fractions of the plankton and the linkage in 
transferring energy from the microbial loop to higher trophic levels 
in the sea [1,2]. Ciliates range in size from <10 to several hundred 
micrometres and may have higher metabolic rates and growth 
efficiencies compared with most metazoan organisms [3]. They 
have a substantial trophic impact in planktonic food webs because 
ciliates can graze 4% to 60% of net annual production, and they 
contribute significantly to heterotrophic biomass and production 
in the aquatic ecosystems. By virtue of their rapid growth rate and 
sensitive reaction to environmental changes, they have been 
considered as effective bioindicators of water quality and 
environmental contamination [2,4,5,6]. 

However, there is still a shortage of data on marine ciliate 
ecology in tropical and subtropical seas up to now [7,8]. Despite 
various investigations on coastal zone resources, ciliate protozoa as 
an important component were usually ignored for their small size. 



fast migration rate, fragile external membranes and the difficulty 
in laboratory culture. Chinese researches on planktonic ciliates 
were mostly conducted in the Yellow Sea [2,4,9,10, 1 1] and South 
China Sea [7,8,12,13,14,15]. In previous studies, ciliate samples 
for enumeration were mostly obtained by net method, which is 
unable to reveal the vertical distribution of ciliates and may cause 
an underestimate in ciliate population [3] . In this study, therefore, 
sampling bottles were used to obtain a more reliable estimate of 
species composition, abundance and biomass of planktonic ciliates. 

The Beibu Gulf is a semi-enclosed gulf located in the northwest 
of the South China Sea surrounded by the land territories of 
China and Vietnam, with the port cities of Beihai, Qinzhou, 
Fangchenggang on the north coast, and Leizhou Peninsula on the 
east. As the fourth largest fish farming area in China, the Beibu 
Gulf has abundant marine resources due to the rich terrigenous 
nutrients carried by coastal rivers. Influenced by the East Asian 
Monsoon, southwesterly winds prevail in summer and northeast- 
erly winds in winter [16,17]. Controlled by multiple water masses 
and currents, the whole Beibu Gulf has complex hydrologic 
condition [17]. There have been a number of investigations on 
plankton community in the Beibu Gulf, however, data on ciliate 
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communities are poor. The spatial pattern of ciliate communities 
has not yet been investigated in this area. We were particularly 
interested in the characteristics of ciliate communities formed in 

specific area. We wonder that of the many environmental factors, 
which are the most tTucial to ciliate communities, although ciliate 
community structure is generally considered as a result of by 
various biotic and abiotic factors. 

The aim of this study was (1) to characterize the species 
composition, abundance, biomass and diversity of the cUiate 
communities from the northern Beibu Gulf, (2) to reveal the 
distribution pattern of the cihate communities at both horizontal 
and vertical levels, and (3) to determine the primary environmental 
factors that controlling the ciliate communities in the northern 
Beibu Gulf in summer. 

Materials and Methods 

Study stations and sampling 

This study was carried out in the northern Beibu Gulf 
(108.25°- 109.88°E, 20.05°~21.46°N), which was surrounded 
by Vietnam, Guangxi and Leizhou Peninsula of China (Figure 1). 
The crui.se was carried out on the vessel "Tianying" from August 
8th to 12th, 2011, and 21 stations were selected in the northern 
Beibu Gulf (Figure 1). Our sampling locations were not national 
park or other protected area of land and no specific permissions 



were required for our locations/activities. We confirm that the 
field studies did not involve endangered or protected species. 

Multilayer water samples were collected from 21 stations using a 
Niskin bottle on the conductivity-temperature-depth (CTD) 
oceanic profilers (SBE-917). At the .stations with water depths 
more than 30 m, water samples were colkx ted at the surface, 
10 m, 30 m, and the bottom, and with water depths less than 
30 m, water samples were collected at the surface, 10 m and the 
bottom, while at the stations with water depths less than 10 m, 
water samples were taken at the surface and the bottom. A total of 
62 water samples were obtained and fixed with acid Lugol's iodine 
solution (final concentration: 2%) in plastic botties. A list of 
sampling layers in 2 1 stations is given in Table 1 . 

Sampling depth, water temperature and salinity were measured 
in situ using an onboard CTD (SBE-917) via probes. Chlorophyll- 
a concentration was measured using a Turner Fluorometer (10- 
AU-005). PH, dissolved oxygen (DO), total organic carbon (TOC), 
nitrate nitrogen (NO3-N), nitrite nitrogen (NO2-N), ammonium 
nitrogen (NH4-N), total nitrogen (TN), active phosphorus (AP), 
total phosphorus (TP) and active silicon (ASi) were determined 
according to the standard methods established in the Offshore 
Marine Chemical Survey Technical Regulations [18]. 

Identification and enumeration 

For identification and enumeration, 1 L of Lugol's fixed 
seawater was settled for 48 hours resulting in 15 mL of 
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Figure 1. Sampling stations of pianittonic ciliates in the northern Beibu Gulf in August 2011. (TIFF). 
doi:l 0.1 371/journal.pone.01 01 121 .gOOl 
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Table 1. Sampling layers in 21 stations for planktonic ciliates 
in the northern Beibu Gulf. 



station 


Sampling layers 


1 surface, 10 m, bottom 


2 


surface, 10 m, bottom 


3 


surface, 10 m, bottom 


4 


surface, 10 m, 30 m, bottom 


5 


surface, 10 m, 30 m, bottom 


6 surface, 10 m, 30 m, bottom 


7 


surface, 10 m, 30 m 


8 


surface, bottom 


9 surface, bottom 


10 


surface, 10 m, bottom 


11 


surface, 10 m, bottom 


12 


surface, 10 m, bottom 


13 


surface, 10 m, bottom 


14 


surface, 10 m, bottom 


15 


surface, 10 m, bottom 


16 


surface, bottom 


17 


surface, bottom 


18 


surface, bottom 


19 


surface, 10 m, bottom 


20 


surface, 10 m, bottom 


21 


surface, 10 m, 30 m, bottom 
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concentrated sample [19]. 0.1 mL of well-mixed concentrated 
sample was taken and identified in a microscope slide. A total 
volume of 1 mL concentrated sample were counted and identified 
under a light microscope at a magnification 400 x . Tintinnids were 
identified by lorica morphology and species description according 
to Hada [20,21,22], Kofoid and Campbell [23,24], Lynn [25] and 
Nie [26]. Aloricate cihateswere identified following Leegaard [27], 
Lynn [25] and Song et al [28]. AH ciliates were finally identified to 
the lowest possible taxa. The individuals which could not be 
identified with Lugol's-fiKed samples were picked out and 
identified using protargol impregnation after re-fixating with 
Bouin's solution [29]. The classification system was mainly 
referred to Lynn [25]. 

Body volumes of ciliate cells were determined from measuring 
of their body size and applying volume equations of approximate 
geometric shape. Conversion factor of C: Vol was 0.19 pg C-|J,m"'' 
for cells preserved with 2% Lugol's iodine [30]. For tintinnids, the 
cell volumes were approximately equal to 1/3 of the loricate 
volumes [8], thus conversion factor of C: Vol was 0.063 pg 
C-|J.m~'' for loricate cells. 

Data analyses 

Species number, abundance and biomass were calculated for 
each sample. The horizontal and vertical distribution patterns of 
planktonic ciliates were displayed using software surfer 8 and 
ocean data view (ODV). 

Shannon-Weiner diversity index [H') for each sample and 
dominance {T) for ciliate species were calculated by the following 
equations: 



s 

H'=-^(ni/N)log2(n,/N) 

i=i 

T= {n,/JV)-f„ Where n, = abundance of the rth species; S = total 
number of species; jV = total number of individuals; f, = 
occurrence of the ith species. 

Depth-averaged abundance and biomass were calculated by the 
integration formula: 

Pd - y — 2^ ^ (A + 1 -JJi) 

Where pj = depth-averaged abundance/biomass; Z^^^^ = 
maximum depth of the station; n = maximum number of layers; 
l<z<jz-l; p, = density in the rth water layer; D, = depth of the rth 
water layer. 

Cluster analysis was used to investigate ciliate community 
structure using the PRIMER v5.0 (Plymouth Routines In 
Multivariate Ecological Research) package according to species 
composition and abundance [31]. To ensure rehabUity, samples 
with observed individuals less than 100 were excluded from the 
biotic data before clustering. Bray-Curtis similarity matrix was 
computed on square root-transformed biotic data and Euclidean 
distance matrix was on ln(x+ 1 )-transformed abiotic data. Group- 
average clustering was used to separate the ciliate communities 
into several assemblages, which represented natural groupings of 
samples (Group A, B, C, etc.), each comprising samples with 
similar species and abundances. Grouping at horizontal level and 
vertical level were also conducted. Groups at horizontal level 
(Group 1, 2, 3, etc.) were classified by clustering based on species 
composition and their depth-averaged abundances in 21 stations, 
and groups at vertical level (Group I, II, III, etc.) were classified 
based on water depth (Group I: 0<water depth<0.5 m. Group II: 
0.5<water depth<10 m. Group III: 10<water depth<30 m. 
Group IV: water depth>30). Variances between groups at 
different levels were tested using the submodule ANOSIM 
(analysis of similarities) in PRIMER v5 (global R values range 
from 0 to 1, and higher global R values are more significant.). 
Species with a cumulative contribution of 90% to the average 
Bray-Curtis similarity within each group was analyzed using the 
submodule SIMPER (Similarity Percentage Analysis). Additional- 
ly, the significance of biota-environment correlation was checked 
by using the RELATE. A potential relationship between biotic 
data and environmental data was explored using the BIOENV 
(biota-environment). To obtain a more detailed relationship 
between ciliate communities and environmental factors, redun- 
dancy analysis (RDA) was conducted by using Caiioco for 
Windows 4.5 package [32]. Under forward selection, significance 
of correlation with environmental variables was determined by 
Monte Carlo permutation test (999 permutations). RDA ordina- 
tion plots were generated based on square root-transformed 
species-abundance data and ln(x+l)- transformed environmental 
data. Spearman correlation analysis between biotic data and 
environmental variables was conducted by using statistical 
program SPSS vl7.0. 
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Results 

Species composition 

A total of 101 species belonging to 44 genera and 7 orders 
(Oligotrichida, Haptorida, Euplotida, Sessilida, Pleurostomatida, 
Scuticociliatida and Tintinnida) were identified. A list of all 
observed species from 7 orders is given in Table SI. The most 
diverse orders were Tintinnida (43 species) and Oligotrichida (40 
species), and the most abundant orders were Oligotrichida 
(2,542.44±248.24 ind-L"') and Haptorida (976.33± 120.82 
ind-L '). As regards ciliate biomass, Tintinnida (2.46±0.65 |J,g 
C-L~') and Oligotrichida (2.31 ±0.25 (Xg C-L"') were the most 
contributive orders. The numbers of loricate and aloricate species 
in all samples were 43 and 58 respectively. Their average 
abundances were 571.81±133.61 ind-L"' and 3,634.96±461.64 
ind-L respectively, and average biomasses were 2.46±0.65 [ig 
C-L"' and 2.73±0.29 ng C-L"', respectively. 

The ciliate communities were dominated by 8 species (domi- 
nance T^0.02), i.e. Mesodinium rubrum (Y-0.\b&2), Strombidium 
inclinatum (3^=0.1018), Mesodinium velox (3^=0.0502), Strombidium 
parastjilifer (2^=0.0481), Strombidium conicum (3^=0.0411), Tintinni- 
dium primitivum {Y= 0.0405), Leegaardiella ovalis {Y= 0.0255), Strombi- 
dium emergens [Y= 0.0229). They occurred in majorit)' of the 
stations, and the cumulative abundance of them exceeded 60% of 
the total in August. Among these dominant species, Mesodinium 
rubrum occurred most frequently that they were observed in 



93.55% of the samples; Tintinnidium primitivum was the only loricate 
species and had a widespread distribution as well. 

Horizontal distribution 

Horizontal distribution of ciliate abundance. In August, 
the ciliate abundance in the surface water in the nearshore area 
was higher than that in the offshore (Figure 2A). The two 
concentrated zones were the coastal region of Guangxi and the 

west coastal region of Leizhou Peninsula, especially the area 
adjacent to Tieshan Harbor. The highest abundance (10.64x10 
ind-L"') was observed in the bottom at station 17. With the 
increasing distance from the coast, the abundance of ciliates 
gradually decreased. The lowest value (0.99x10^ ind-L"') was 
observed at station 3. The distribution in the bottom layer showed 
a similar trend with surface water (Figure 3A). The highest 
abundance (10.95 x 10 ind-L ) occurred in the bottom at station 
18, while the lowest value (0.74x10'' ind-L ') occurred at station 
3. Concerning the depth-averaged abundance, the high values was 
observed around stations 17 and 18, while the lowest value was 
recorded at station 3 (Figure 4A). 

Horizontal distribution of ciliate biomass. The depth- 
averaged biomass decreased westwards from the west coast of 
Leizhou Peninsula (Figure 4B). The same was true for the biomass 
in each layer. In the surface water, the highest biomass (18.26 \xg 
C-L ) was observed at station 18 (Figure 2B), while the lowest 




Figure 2. Horizontal distribution of planktonic ciliates in the surface layer of the northern Beibu Gulf. (Aj Horizontal distribution of 
abundance. (B) Horizontal distribution of biomass. (C) Horizontal distribution of species number. (D) Horizontal distribution of diversity. The 
horizontal axis and the vertical axis expressed the gradients in longitude and latitude, respectively. (TIFF). 
doi:l 0.1 371/journal.pone.01 01 121 .g002 
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Figure 3. Horizontal distribution of planlctonic ciiiates in the bottom layer of the northern Beibu Gulf. (A) Horizontal distribution of 
abundance. (B) Horizontal distribution of biomass. (C) Horizontal distribution of species number. (D) Horizontal distribution of diversity. The 
horizontal axis and the vertical axis expressed the gradients in longitude and latitude, respectively. (TIFF). 
doi:l 0.1 371/journal.pone.01 01 121 .g003 



value (0.48 \lg C-L ') occurred at station 3. In the bottom water, 
the ciliate biomas.s pealced at station 18 (30.53 \xg C-L '), and the 
lowest value (0.45 ng C-L ') occurred at station 1 (Figure SB). 

Horizontal distribution of species number. The species 
numbers of ciiiates were variable among stations. In the surface 
water, a conspicuous westwards decrease of the species numbers 
was showed, ranging from 3 1 (station 1 2) in the east to 8 (station 3) 
in the west (Figure 2C). The bottom layer showed a similar trend 
that the species numbers in the west area were less than in the east 
(Figure 3C). For the species number in each station, however, 
higher values were concentrated in the area north of Hainan 
Lsland while lower values were recorded in the coastal waters of 
Guangxi- Vietnam (Figure 4C). There were only 18 species 
identified at station 1, while as many as 47 species were observed 
at station 19. 

Horizontal distribution of diversity. In contrast with the 
distribution trend of ciliate abundance (Figure 2A, 3A, 4A), 
Shannon- Weiner diversity of each station was obviously lower in 

the nearshore waters than in the offshore waters (Figure 4D), and 
the same was true for the diversity in each layer (Figure 2D and 
3D). The zones with lower values displayed patchy distributions. 
One patch was around station 8 at the mouth of Qinzhou Bay in 
Guangxi, and the other was around station 20 in the west coast of 
Leizhou Peninsula. The area nearby Tieshan Harbor as well as the 
large area southwest of Weizhou Island displayed higher values. In 
the surface water. Shannon- Weiner diversity minimized (1.82) at 



station 1 and maximized (4.17) at station 12 (Figure 2D). In the 
bottom layer, the minimum (1.81) was observed at station 20, 
while the maximum (4.35) was at station 13 (Figure 3D). 

Vertical distribution 

In this study, four vertical sections were picked out to investigate 
the vertical distribution trend of ciliate abundance, biomass, 
species number and Shannon-Weiner diversity with increased 
water depth and distance from the coast. Sections 1^ (stations 1, 
2, 3, 4) and 5-8 (stations 5, 6, 7, 8) were in longitude direction, 
while sections 4-19 (stations 4, 5, 12, 13, 19) and 3-20 (stations 3, 
6, 11, 14, 20) were in latitude direction. 

Along the section 1-4, abundance of ciiiates were low in the 
surface water but it maximized (8.53x10^ ind-L ') at a depth of 
10 m in the nearshore waters (Figure 5A). The second highest 
value was observed in the layer of 30 m in the offshore waters. The 
distribution of the biomass, in accord with abundance, showed 
high values at a depth of 10 m in the nearshore waters and at 
about 30 m offshore (Figure 5B). In the section 5-8, the 
abundance peaked with a higher value 9.92x10'^ ind-L~' at a 
depth of 10 m (Figure 6A). The high values of biomass, however, 
only occurred at a depth of 30 m in the offshore waters (Figure 6B). 
The abundance of ciiiates in the section 4-19 showed a maximum 
value (6.83x10^ ind-L"') at a depth of 30 m in offshore waters, 
and higher values were also observed at the surface layer in the 
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Figure 4. Horizontal distribution of planlctonic ciliates in tlie water columns of the northern Beibu Gulf. (A) Horizontal distribution of 
abundance. (B) Horizontal distribution of biomass. (C) Horizontal distribution of species number. (D) Horizontal distribution of diversity. The 
horizontal axis and the vertical axis expressed the gradients in longitude and latitude, respectively. (TIFF). 
doi:l 0.1 371/journal.pone.01 01 121 .g004 
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Figure 5. Vertical distribution of planlttonic ciliates in the section 1-4. (A) Vertical distribution of abundance. (B) Vertical distribution of 
biomass. (C) Vertical distribution of species number. (D) Vertical distribution of diversity. The horizontal axis and vertical axis expressed the gradients 
in latitude and water depth, respectively. (TIFF). 
doi:1 0.1 371/journal.pone.01 01 121 .g005 
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Figure 6. Vertical distribution of planktonic ciliates in the section 5-8. (A) Vertical distribution of abundance. (B) Vertical distribution of 
biomass. (C) Vertical distribution of species number. (D) Vertical distribution of diversity. The horizontal axis and vertical axis expressed the gradients 
in latitude and water depth, respectively. (TIFF). 
doi:1 0.1 371/journal.pone.01 01 121 .g006 



nearshore waters (Figure 7A). The biomass maximized at a depth 
of 10 m in the nearshore waters and the second highest biomass 
occurred at a depth of 30 m in the offshore waters (Figure 7B). In 
the section 3-20, the highest abundance and biomass were both 
occurred at a depth of about 30 m in offshore waters. 

High species numbers in the section 1-4 gathered at depths of 
10 m and 30 m in the offshore area (Figure 5C). In the section 5— 
8, high species numbers appeared below the depth of 10 m in the 
offshore waters (Figure 6C). In the section 4—19, the higher values 



Abundance pnd/L] 




20.4'H 20.rN 20.2>M 20.1'N 



Species number 




20.4'N 20.rN 20.2'N 20.1'N 



were mostly observed in the nearshore waters, and the values were 
also relatively high at a depth of more than 20 m in the offshore 
waters (Figure 7C). In the section 3-20, lower species numbers 
were recorded in the nearshore waters while higher values were 
observed at about 10 m in the offshore waters (Figure 8C). Areas 
with higher Shannon- Weiner diversity indices were almost those 
with higher species numbers (Figure 5-8). 

A conspicuous consistent trend could be found out between the 
vertical distributions of ciliate abundance and biomass. The 



Biomass [/ig C/L] 




20.4'N 20.3*M 20.2*M 20.1'N 



H' 
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Figure 7. Vertical distribution of planktonic ciliates in the section 4-19. (A) Vertical distribution of abundance. (B) Vertical distribution of 
biomass. (C) Vertical distribution of species number. (D) Vertical distribution of diversity. The horizontal axis and vertical axis expressed the gradients 
in latitude and water depth, respectively. (TIFF). 
doi:1 0.1 371/journal.pone.01 01 121 .g007 
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Figure 8. Vertical distribution of planktonic ciliates in the section 3-20. (A) Vertical distribution of abundance. (B) Vertical distribution of 
biomass. (C) Vertical distribution of species number. (D) Vertical distribution of diversity. The horizontal axis and vertical axis expressed the gradients 
in latitude and water depth, respectively. (TIFF). 
doi:1 0.1 371/journal.pone.01 01 121 .g008 



aggregated areas were displayed at a depth of 10m in the 

near.shore waters, while descended to about 30 m in the ofTshore 
waters. No consistent trend was found between species number 
and abundance/biomass, however, the vertical distributions of 
species number and diversity were almost the same. 

Spatial pattern of ciliate communities 

The variations of horizontal distribution and vertical distribu- 
tion were studied. At the horizontal level, cluster analysis on the 
data of species composition and their depth-averaged abundance 
divided 2 1 stations into four groups at a similarity level of 50. 1 1 % : 
Group 1 comprised 3 stations (18, 20, 21) located in the west coast 
of Leizhou Peninsula; Group 2 included only one station (17) 
situated nearby Tieshan Harbour; Group 3 covered most oflF-shore 
stations (2-7, 9-16, 19) in the southwest of our study area; Group 
4 included 2 inshore stations (1, 8) by the coast of Guangxi. Based 
on the groupings above, all 62 samples were classified into four 
groups. A relatively significant difference between groups was 
tested at horizontal level (global i? = 0.376, P<0.001). At the 
vertical level, all samples were also divided into four groups in 
terms of the depth of each sample. Group I consisted of 21 samples 
from a depth not more than 0.5 m; Group II comprised 20 
samples from a depth more than 0.5 m but not more than 10 m; 
Group III comprised 17 samples from a depth more than 10 m 
but not more than 30 m, and Group IV included the other 4 
samples which were collected from a depth more than 30 m. The 
difference between groups was not so significant (global R — 0.095, 
P= 0.002), suggesting a less conspicuous variation of the vertical 
distribution than of the horizontal distribution. 

Cluster analysis separated the ciliate samples in which 
individuals were not less than 100 into several groups, each 
comprising samples with similar species composition and biomass. 
A dendrogram showed that 5 1 effective samples were divided into 
four clusters at a similarity level of 27.51% (Figure 9). Thus three 
main groups were obtained excluding the bottom sample of station 
2. Group A was comprised of samples of stations 18-21 west of 



Leizhou Peninsula and at the mouth of Qiongzhou Strait; Group 

B consisted of off-shore stations 4, 5, 6, 1 1 in the southwest area; 
Group C consisted of the most samples from the southwest off- 
shore stations. Group C was further divided into three subgroups. 
Subgroup a consisted of samples of stations 15, 17 nearby Tieshan 
Harbour and of station 2. Subgroup b consisted of samples of 
stations 1 and 8. Subgroup c consisted of the other samples from 
the southwest area. 

Species that were most contributive to each group were 
determined using SIMPER analysis. In Group A, Tintimidium 
primitimm, the only tintinnid of the dominant species, was with a 
greatest contribution of 44.28%, followed by Ttntinnopsis schotti 
(19.34%) and Laboea strobili (13.78%). Major contributors to Group 
B were Leegaardiella ovalis (20.72%), Tintinnidium primitimm (17.00%) 
and Strombidium pamcalkinsi (16.81 '^i)). Group C, comprising of most 
samples, was the largest group, and the most contributive species 
to Group C were Omegastrombidium jankowskii (14.00%), Strombidium 
conicum (11.67%), Mesodinium rubrum (10.82%). Within Group C, 
Omegastrombidium jankowskii (32.34%) and Strombidium conicum 
(10.57%) were most contributive to Subgroup a, Strombidium 
inclinatum (51.71%) and Strombidium conicum (15.93%) were most 
contributive to Subgroup b, whilst Omegastrombidium jankowskii 
(12.78%) and Mesodinium rubrum (10.73%) were more important in 
Subgroup c. The bottom sample from station 2 was not included 
in any of these groups, mainly because of the dominance of 
Leprotintinnus simplex, Stenosemella pawicollis and Strombidium conicum. 

Correlations between ciliate communities and 
environmental factors 

RELATE analysis revealed a significant correlation between the 
spatial variation of ciliate communities and the variation of 
environmental data (fl = 0.109; P= 0.004). BIOENV analysis 
suggested the combination of environmental variables that were 
most correlative with the spatial variation of ciliate communities 
were: AP, TP, NOg-N, NO3-N and NH4-N (Table 2). Among 
them TP was the only variable included in all the top 10 matches, 
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indicating TP was the key environmental variable that significantly 
influenced the variation of ciliate communities. Spearman 
correlation analysis indicated that the ciliate abundance, biomass 
and species numbers were significantly and positively correlated 
with concentration of chlorophyll-a, and ciliate biomass also 
exhibited a significantly positive correlation with concentrations of 
phosphorus and nitrogen nutrients (AP, TP, NO2-N) (Table .3). 

A species-abundance matrix in detrended correspondence 
analysis (DCA) resulted in the biggest gradient length 3.320. 
Redundancy analysis (RDA) method was selected for the analysis 



of sample-environment relationship. Among 14 environmental 
factors for study, four most significant factors, i.e. TP, NO3-N, 
chlorophyll-fl and water temperature, were picked out by forward 
selection. All the four canonical axes jointly explained 77.7% 
(significant factors) of the species-environment relationship. The 
first and second axes which jointly explained 63.5% of the 
relationship were exhibited on the RDA two-dimensional plot 
(Figure 10). Group A was at a considerably high level of TP and 
NO3-N. Group B showed lower levels of water temperature and 
chlorophyll-a. By contrast, Group C displayed a higher level of 
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Table 2. Result of biota-environment (BIOENV) analysis. 



Rank 


R 


Environmental variables 


P 


1 


0.435 


AP, TP, NO2-N 


0.01 


2 


0.427 


TP, NO2-N 


0.01 


3 


0.407 


AP, TP, NO2-N, NH4-N 


0.01 


4 


0.404 


AP, TP, NO2-N, NO3-N 


0.01 


5 


0.400 


TP, NO2-N, NH4-N 


0.01 


6 


0.399 


TP, NO2-N, NO3-N 


0.01 


7 


0.391 


AP, TP, NO2-N, NO3-N, NH4-N 


0.01 


8 


0.387 


TP, NO2-N, NO3-N, NH4-N 


0.01 


9 


0.383 


AP, TP 


0.01 


10 


0.380 


AP, TP, NO3-N 


0.01 



The result shows the top 10 matches of environmental variables with spatial 

variance of ciliate community. 

doi:l 0.1 371 /journal.pone.Ol 01 1 21 .t002 



water temperature and chlorophyll-a. In Group C, there was 
difference between subgroups. Water temperature and chloro- 
phyll-fl concentration in Subgroup a and c were relatively high, 
whilst those in subgroup b were lower. 

Discussion 

Spatial pattern of ciliate communities 

A high value of diversity index suggests a more healthy 
ecosystem, while a low value suggests a less healthy or degraded 
ecosystem [33]. Shannon- Weiner diversity index H' is one of the 
most widely used diversity indices for measuring biodiversity. The 
Shannon-Weiner diversity index of ciliate communities in the 
northern Beibu Gulf remarkably varied. It was much lower in 
coastal waters than in offshore waters, implying a great impact of 
human activity to the water quality and the health of plankton 
ecosystem. 

In addition, ANOSIM analysis revealed that the variation of 
cUiate communities was significant at horizontal level, but it was 
not significant at vertical level. The strong tide in the Gulf gives the 
basic force for vertical and horizontal mixing, making the waters 
mix quickly [16]. At vertical level, the water is relatively shallow 
that the maximum depth is less than 100 m. Seawater is well 
mixed in the vertical direction so that the physical and chemical 
conditions in the whole water columns vary slightly. At horizontal 
level, in contrast, the northern Beibu Gulf covers a vast expanse 
and is affected by different water masses and currents. The 
physical and chemical characters of seawater and how badly it is 
affected by human activities differ gready. As a result, the ciliate 
communities in the study area exhibited a relatively homogeneous 
distribution vertically but varied gready among stations. 

Relationship between ciliate communities and 
environmental factors 

Tintimidium primitivum is a typical species living in the cold water 
mass under the thermocline. In our survey, T. primitivum did 
existed in the cold water mass in the southwest area, verifying its 
positive correlation with water depth and negative correlation with 
temperature (Table 3). However, a much greater abundance of T. 
primitivum was recorded in the coastal waters west of Leizhou 
Peninsula. T. primitivum was also proved to be positively correlated 
with nitrogen and phosphorus (AP, TP, NOg-N, NO3-N) (Table 3), 
revealing its preference for nitrogen and phosphorus. 



Among the dominant species, two [Mesodinium ruhmm and 
Mesodinium velox) were in genus Mesodinium and they contributed 
significandy to the total abundance. Mesodinium spp. can form a 
high density in a short time in eutrophic waters and cause 
Mesodinium spp. bloom [7]. M. rubrum is an autotrophic species 
containing a permanent symbiont consisting of many chloroplasts 
[34-36]. It can utilize light and nutrients to carry out photosyn- 
thesis and has strong ability of motility and vertical migration [35]. 
The highest abundance (4.46x10'^ iiid-L ') oi M. rubrum was 
observed in the surface waters west of Leizhou Peninsula, where a 
high concentration of chlorophyU-a was measured. Our results 
showed that M. rubrum was significantly correlated with total 
chlorophyll-fl as well as nano- and micro- sized chlorophyU-a 
concentration. The contribution of M. rubrum to the chlorophyU-a 
concentration may be important. And this may also suggest a 
preference oi M. rubrum to eutrophic conditions. 

The abundance, biomass and species number of ciliate 
communities were found significantly correlated to chlorophyll-a 
concentration, especially the nano- and micro- chlorophyU-a. For 
one thing, heterotrophic planktonic ciliates respond rapidly to 
their food availability in the aquatic environment. For another, 
autotrophic ciliates (e.g. Mesodinium rubrum) and mixotrophic 
ciliates (e.g. Laboea strobila, Strombidium stylifef) contain chloroplasts 
or symbiont consisting of chloroplasts [34—38]. They may 
contribute greatly to the chlorophyU-a in the northern Beibu 
Gulf, as in many coastal and estuarine ecosystems [39-41]. 

Besides, the spatial variation of ciliate communities was 
significandy correlated with the variables of phosphorus and 
nitrogen nutrients (i.e. TP, AP, NOj-N and NO3-N), suggesting a 
potential impact of eutrophication in estuarine and coastal waters 
on the ciliate communities. Nutrients direcdy affect the photosyn- 
thesis of photosynthetic ciliates, while heterotrophic ciliates are 
impacted by nutrients indirecdy, through the effect on phyto- 
plankton growth. Spearman analysis also indicated high concen- 
trations of phosphorus and nitrogen nutrients primarily influenced 
ciliate biomass. It means the eutrophic state is more likely to cause 
an increase in ciliate size rather than their quantity or diversity. 

The ciliate community is affected by multiple environmental 
factors. Both top-down and bottom-up interactions may play a 
role in population control of ciliates [42] . The response of ciliate 
community to nutrients and chlorophyU-a concentration also 
provided an evidence for ciliates as bioindicators. 

Relationship between ciliate communities and water 
masses 

It should be noted that planktonic ciliates are creatures that lack 
of ability of independent movement and drift with the water 
currents. Their spatial distribution largely depends on the water 
currents and masses in study area. The whole Beibu Gulf can be 
generally divided into several masses, i.e. the coastal water, the 
outer Beibu Gulf water, the mixed water, the bottom water and so 
on [16,17]. The coastal water mass with low salinity mainly affects 
the northwest of the bay, and the South China Sea water mass 
with high salinity covers the south of the bay. Meanwhile, the 
northern Beibu Gulf is also influenced by coastal currents and 
Qiongzhou Strait West Current. 

In the Beibu Gulf, the southwesterly wind strengthens the 
northward flow west of Hainan Island in summer, causing the 
strongest intrusion of the South China Sea surface water from the 
mouth of the gulf, which is characterized by high salinity [17]. 
Affected by the surface and subsurface water, the bottom water 
which originates from the coastal water is limited below the 
thermocline and covers only a small area east of BaUongwei Island 
[16,43]. In our study, the thermocline was observed at a depth of 
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l(}-30 m in stations 4, 5, 6 and 11. Group B may represent an 
ofTshore community which was mainly controlled by the bottom 
cold water. Leegaardiella ovalis, a small sized species, which prefers 
high salinity and low temperature, prevailed in this group. 

The diluted water, characterized by low salinity and relatively 
high temperature and dissolved oxygen, is mainly formed by the 
river discharges from the Nanliu River, the BeUun River, and the 
Dafeng River in Guangxi. Interacted with the high-salinity South 
China Sea water, the diluted water is restricted to a small area off 
the Guangxi- Vietnam coast. Most area north of 20°N and west of 
Leizhou Peninsula was occupied by the mixed water [16]. The 
mixed water is a mixture between the diluted water and the 
surface water with salinity and temperature intermediate between 
them. The chlorophyll-a concentration in the mixed water is 
apparently higher than that in the surface water or the subsurface 
water. Group C covered quite a large area in the west of our study 
area. It was influenced by not only the diluted water but also the 
surface water, showing a preference for high water temperature 
and chlorophyll-fl concentration. Larger naked species, Omegas- 
trombidiiim jankowskii and Strombidium conicum, were dominant in 
Group C. Meanwhile, Mesodinium rubrum, as a typical autotroph, 
was abundant in this group. The bottom water in station 2, 
different from other stations nearby, was characterized by low 
levels of temperature and chlorophyll-a. The species composition 
here was distinctive. The special nature of station 2 was speculated 



to be resulted from the diluted water from Vietnam, however, 
more evidence was needed for further confirmation. 

The outer surface water is pushed northward along the western 
coast of Hainan Island after coming into the gulf [16]. It joins the 
Qiongzhou Strait West Current and then flows northward along 
the western coast of Leizhou Peninsula. A high level of nutrients, 
which may be brought by Qiongzhou Strait West Current from 
the east coast of Leizhou Peninsula, was detected in the waters 
west of Leizhou Peninsula. Group A in this area showed a 
preferrence to high concentration of nutrients and was dominated 
by loricate species Tintinnidium primitivum. 

Conclusion 

The ciliate communities in summer in the northern Beibu Gulf, 
South China Sea showed significant spatial variations of species 
composition, abundance, biomass and diversity. The variation of 
ciliate communities was significant at horizontal level, but that was 
not at vertical level. The spatial pattern of ciliate communities was 
closely related to certain environmental factors, among which 
phosphorus, nitrogen and chlorophyll-a were the controlling 
factors in summer in this area. Our results will provide basic data 
for further studies in ciliate community structure and in aquatic 
ecology in the Beibu Gulf 



PLOS ONE I www.plosone.org 



12 



July 2014 I Volume 9 | Issue 7 | e101121 



Factors Controlling Ciliate Communities 



Supporting Information 

Table SI List of the ciliate species recorded in multi- 
layer samples &om the northern Beibu Gulf. 

(DOC) 

Aclcnowledgments 

We thank Prof. Henglong Xu from Ocean University of China for his help 
on identifying species and Dr. Lemian Liu from Institute of Urban 

References 

1. Azam F, Fenchel T, Field JG, Gray JS, Meyer-Reil LA, et al. (1983) The 
ecological role of water-column microbes in the sea. Mar Ecol Prog Ser 10: 257- 
263. doi: 10.3354/mcps010257. 

2. Jiang Y, Xu H, Al-Rasheid KAS, Warren A, Hu X, et al. (2011) Planktonic 
ciliaTc communidrs in a semi-enclosed bay of Yellow Sea, northern China: 
annual cNdc. Jo Mar Biol Ass6 U6 K 91: 97-105. doi; 10.1017/ 
S002r)3lr)41000175X. 

3. Montagncs DJS, Lynn DH, Rofi JC, laylor VVU (1988) The annual cycle of 
heterotrophic planktonic ciliates in the waters surrounding the Isles of Shoals, 
Gulf of Maine: an assessment of their trophic role. Mar Biol 99: 121-30. doi: 
10.1007/BF00644973. 

4. Jiang Y, Xu H, Hua X, Zhu M, Al-Rasheid KAS, et al. (201 1) An approach to 
analyzing spatial patterns of planktonic ciliate communities for monitoring water 
quality in Jiaozhou Bay, northern China. Mar PoUut Bull 62: 27-235. doi: 
1 0. 1 0 i G/j.marpolbul.20 1 0. 1 1 .008. 

5. Kim Y()^ Shin K, Jang PC, Choi HW, Noh JH, et al. (20 1 2) I'intinnid species as 
biological indicators for monitoring intrusion of the warm oceanic waters into 
Korean coastal waters. Ocean Sci5 J 47: 161-172. doi: 10.1007/sl2601-012- 
0016-4. 

6. Xu H, Song W, Warren A, Al-Rasheid KAS, Al-Farraj SA, et al. (2008) 
Planktonic protist eommtmities in a semi-enclosed mariculture pond: structural 
variation and correlation with enx ironmental conditions. J Mar Biol Ass UK 88: 
1353-1362. doi: 10.1017/S0025315408002129. 

7. Liu H, Tan Y, Huang L, Song X, Huang J, et al. (2010) Composition and 
distribution of ciliates in northern South China Sea during summer. Acta Leol 
Sin 30: 2340-2346. 

8. Zhang C, Zhang W, Xiao T (2010) Ciliate abundance and biomass in northern 
Soudi China Sea in October 2007. Acta Ecol Sin 30: 0867-0877. 

9. Song W, Wang M, Warren A (2000) Redescriptions of three marine ciliates, 
Strombidium degans Florentin, 1901, Stwmbidium sulcatum Claparede & Lachmann, 
1859 and HeUmstrombidium paracalkinsi Lei, Xu & Song, 1999 (Ciliophora, 
Oligotrichida). Europojo Protistology 36: 327-342. 

10. Xu D, Song \V (2006) Hapantotvpiliralion and morjihological redeseription of 
the marine planktonic ciliate, Spirostrombidium cinduni (Kahl, 1932) Petz, Song et 
Wilbert, 1995 (Ciliophora: Oligotrichida). Acta Protozool 45: 17-25. 

1 1 . Zhang W, Xu K, Wan R, Zhang G, Meng T, et al. (2002) Spatial distribution of 
ciliates, copepod nauplii and eggs, Etigrmdis japonicus post-larvae and micro- 
zooplankton herbivorous activity in the Yellow Sea, China. Aquat Microb Ecol 
27: 249-259. doi: 10.3354/ame027249. 

12. Feng M, Zhang W, Zhang C, Xiao T, Li C (2010) Horizontal distribution of 
large tintinnids in the northern South China Sea. J IVop Oceanogr 29: 141-150. 
doi: 10.3969/J.issn.l009-5470.2010.03.023. 

13. Gomez F (2007) Trends on the distribution of ciliates in the open Pacific Ocean. 
Acta Oecol 32: 188-202. doi: 10.1016/j.actao.2007.04.002. 

14. Su Q, Huang L, Tan Y, Xu R, Li T, et al. (2007) Preliminary study of 
microzooplankton grazing and community composition in the north of South 
China Sea in autumn. Mar Sci Bull 9: 43-53. doi: 10.3969/j.issn.lOOO- 
9620.2007.02.006. 

15. Xu R, Bai 0^(1998) Strombidinopsis grmdis: A new species of Strombidinopsis. Trop 

Oeeanol 17: 40-43. 

16. Chen S, Li Y, Hu J, Zheng A, Huang L, et al. (2011) Multiparameter cluster 
analysis of seasonal variation of water masses in the eastern Beibu Gulf. 
Jo Oceanography 67: 709-718. doi: 10.1007/sl0872-011-0071-y 

17. Sun X (2005) The hydrography in the Beibu Gulf In: SuJ, editor. Hydrography 
in the China coastal ocean. Beijing: Ocean Press, pp. 285. 

18. Office of the State Oceanic Administration (2006) Comprehensive survey and 
assessment technique regulations in offshore area of China. Marine chemistry 
survey technique regulations. Beijing: Ocean Press. 

19. Utermohl H (1958) Zurvervolkommungder quantitativen phytoplankton Meth- 
odik. Mitt Int Verein Theor Angew Limiiol 9: 1-38. 

20. Hada Y (1932) The Tintinnoinea from the sea of Okhotsk and its neighborhood. 
J6 Fac Sci Hokkaido Imp Univ ser VI, Zool 2: 37-60. 



Environment, Chinese Academy of Sciences for comments on an earlier 
draft of the manuscript. We also thank Marine Chemistry Laboratory 
(Prof. Airong Zheng), Xiamen University, China for providing chemical 
data. 

Author Contributions 

Conceived and designed the eijperiments: WZ. Performed the experiments: 
YW WZ. Analyzed the data: YW JY. Contributed reagents/materials/ 
analysis tools: WZ YL WC LZ. Wrote the paper: YW WZ JY. 



21. Hada Y (1937) The fauna of Akkeshi Bay, W. The pelagic Ciliata. Ja Fac Sci 
Hokkaido Imp Univ ser FV, Zool 5: 143-216. 

22. Hada Y (1938) Studies on the Tintinnoinea from the western tropical Pacific. 
Jo Fac Sci Hokkaido Imp Univ ser FV, Zool 6: 87-190. 

23. Kofoid CA, Campbell AS (1929) A conspectus of the marine and freshwater 
Ciliata belonging to the suborder Tintinnoinea, with descriptions of new species 
principally from the Agassiz expedition to the eastern tropical Pacific 1904- 
1905. Univ Calif Publ Zool 34: 1-403. 

24. Kofoid CA, Campbell AS (1939) I'he I'intinnoinea of the eastern tropical 
Pacific. Bull .\lus Comp Zool Harvard Coll, 84: 1-473. 

25. Lynn DH (2008) The ciliated protozoa. Characterization, classification and 
guide to the literature, 3rd eel. Springer, Berlin. 

26. Nie D (1934) Notes on Lintinnoinea from the Bay of Amoy. Ann Rep Mar Biol 
Assoc China 3: 71-80. 

27. Leegaard C (1915) Untersuchungen tiber einige planktonciliaten des meeres. 
Nyt Mag f Naturv 53: 1-37. 

28. Song W, Warren A, Hu X (2009) Free-living ciliates in the Bohai and Yellow 
Seas, China, Beijing: Science Press, pp. 1—518. 

29. Montagnes DJS, Humphrey E (1998) A Description of Occurrence and 
Morphology' of a New Species of Red- Water Forming Strombidium (Spiro- 
trichea, Ohgotrichia). J Euk Microbiol, 45: 502-506. doi: 10.1111/j.l550- 
7408.1998.tb05108.x. 

30. Putt M, Stoecker DK (1989) An experimentally determined carbon: volume 
raUo for marine "oligotrichous" ciliates from estuarine and coastal waters. 
Limnol Oceanogr 34: 1097-1103. doi: 10.4319/lo.l989.34.6.1097. 

31. Clarke KR, Gorley RN (2001) PRIMER v5: Usermanual/tutorial. Plymouth, 
UK; PRIMER-E. 

32. ter Braak CJF, Smilauer P (2002) CANOCO Reference Manual and CanoDraw 
for Windows User's Guide: software for canonical community ordination 
(version 4.5). Microcomputer Power, Ithaca, New York. 

33. Wilhm JL (1970) Range of diversity index in benthic macroinvertebrate 
populations. Jo Water Poll 42: 221-224. 

34. Stoecker DK, Silver MW (1987) Chloroplast Retention by Marine Planktonic 
Ciliatcsa. Ann NY Acad Sci, 503; 562-565. doi: 10.1111/j.l749- 
6632.1987.tb40646.x. 

35. Crawford DW (1989) Mesodinium nibrum: the phytoplankter that wasn't. Mar Eco 
Prog Ser 58: 161-174. doi: 10.3354/meps058161. 

36. Gustafson DE, Stoecker DK, Johnson MD, Van Hetikelem WF, Sneider K 
(2000) Cryptophyte algae are robbed of their organelles by the marine ciliate 
Mesodinium mbmm. Nature 405: 1049-1052. doi; 10.1038/35016570. 

37. Dolan JR, P^rez MT (2000) Costs, benefits and characteristics of mixotrophy in 
marine oligotrichs. Freshw Biol 45: 227-238. doi: 10.1046/j.l365- 
2427.2000.ob659.x. 

38. McManus (iB, Zhang H, Lin SJ (2004) Marine planktonic cihates that prey on 
maeroalgae and enslave their chloroplasts. Limnol Oceanogr 49: 308—313. doi: 
10.4319/10.2001.19.1.0308. 

39. Smith \\(). Barber RT (1979) Carbon budget for the autotrophic ciliate 
Mesodimimimbmin.} Phycol 15: 27-33. doi: 10.1 1 1 l/j.l529-8817.1979.tb02957.x. 

40. Stoecker DK, Putt M, Da\is LH, Michaels AE (1991) Photosynthesis in 
Mesodinium rubnm: species-specific measurements and comparison to community 
rates. Mar Ecol Prog Ser 73: 245-252. doi:10.3354/meps073245. 

41. Johnson MD, Stoecker DK, Marshall HG (2013) Seasonal dynamics of 
Mesodinium rubrum in Chesapeake Bay. Jo Plankton Res 35: 877—893. 
dot: 10. 1 093/plankt/flit028. 

42. Sanders RW (1995) Seasonal distributions of the photosvnthesizing cihates 
Lfihoea slwbila and Mi rionfcta rubra ( — Mesodinium rubrum) in an estuary of the Gulf 
of Maine. Aquat Microb Ecol 9: 237-242. doi: 10.3354/ame009237. 

43. Chen S, HuJ, ZhtiJ, Sun Z, Chen Z (2008) Water mass analysis in the eastern 
Beibu Gulf in summer. In: HuJ, Yang S, editors. A collection of research papers 
on marine science in the Beibu Gulf, vol 1. Ocean, Beijing, pp. 88-97. 



PLCS ONE I www.plosone.org 



13 



July 2014 I Volume 9 | Issue 7 | e101121 



